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Abstract. We consider non-negative solutions of the fast diffusion equation u t = Au m with 
m 6 (0, 1), in the Euclidean space R d , d > 3, and study the asymptotic behavior of a natural 
class of solutions, in the limit corresponding to t — > oo for m > m c = (d — 2)/d, or as t approaches 
5— I ■ the extinction time when m < m c . For a class of initial data we prove that the solution converges 

with a polynomial rate to a self-similar solution, for t large enough if m > m c , or close enough 
to the extinction time if m < m c . Such results are new in the range m < m c where previous 



00 

< 



approaches fail. In the range m c < m < 1 we improve on known results. 



1. Introduction 



We study the Cauchy problem for the fast diffusion equation posed in the whole Euclidean 



space, that is, we consider the solutions u(t, y) of 

( d T u 



1.1 



u(0, •) = u 



> . 

CN , where m G (0, 1) (which means fast diffusion) and (r, y) E (0, T) x IT for some T > 0. We 

consider non-negative initial data and solutions. Existence and uniqueness of weak solutions 
of this problem with initial data in Lli oc (lR' i ) was first proved by M.A. Herrero and M. Pierre 
in [30] . In the whole space, the behavior of the solutions is quite different in the parameter ranges 
■ m c < m < 1 and < m < m c , the critical exponent being defined as 

' d-2 
O ■ m c := — — . 

a 

Note that m c > only if d > 3, so that the lower range does not exist for d = 1, 2. For m > m c 
^ ' the mass J Rd u(y, t) dy is preserved in time if the initial datum uq is integrable in M. d . Besides, 

non-negative solutions are positive and smooth for all x £ M. d and t > 0. On the contrary, 
solutions may extinguish in finite time in the lower range m < m c , for instance when the initial 
data is in L p *(]R d ) with = d (1 — m)/2: then there exists a time T > such that 

lim u(t, y) = . 

Many computations are however similar in both ranges, from an algebraic point of view. We 
refer to the monograph [47] for a detailed discussion of the existence theory and references to the 
subject. The extension to exponents m < is also treated, and it is natural but it will not be 
the focus of this paper. 

In the last two decades, special attention has been given to the study of large time asymptotics 
of these equations, starting with the pioneering work of A. Friedman and S. Kamin [28] and 
completed in [35], when m is in the range (m c ,oo). In those studies the class of non-negative, 
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finite mass solutions are considered. Asymptotic stabilization towards self-similar asymptotic 
solutions known as Barenblatt solutions is shown. For m c < m < 1, such solutions take the form: 



;i-2) U D , T (T,y) :-- 




R(t) 



1 

1 — m 



with R{t) := [d ( m — m c ) (r + T)] d ( m -™c) . Here D, T > are free parameters. While the second 
parameter means a time displacement and does not play much role in the asymptotic behavior, 
the first does and can be computed from the mass of the solution. The value m c is the critical 
exponent below which the Barenblatt solutions cease to exist in this standard form. 

Here, we are mainly interested in addressing the question of the asymptotic behavior of 
when < m < m c . We consider a wide class of solutions which vanish in finite time T and describe 
their behavior as r goes to T. We point out that our methods allow to treat simultaneously the 
ranges < m < m c and m c < m < 1, in which one is interested in the behavior of the 
solutions as r goes to infinity. For this purpose, we extend the Barenblatt solutions to the range 
< m < m c with the same expression (jl.2p . but a different form for R, that is 



R(t) := [d (m c — m) (T — r)] ^^7. 

The parameter T now denotes the extinction time. Following [U], we shall call such solutions 
the pseudo-Barenblatt solutions. Notice that Barenblatt and pseudo-Barenblatt solutions Ud,t, 
with D, T > 0, are such that \J V D T is integrable if and only if p > p* (p* is defined above, and 
p* > 1 means m < m c ). Consistently with the above choices, for m = m c , one has to choose 
R(t) := e T+T with free parameter T, see |47| . in order to obtain pseudo-Barenblatt solutions; 
then, p* = 1. 

The family of Barenblatt (respectively pseudo-Barenblatt) solutions represents the asymptotic 
patterns to which many other solutions converge for large times if m > m c (respectively as t goes 
to T if < m < m c ). We are interested in the class of solutions for which such a convergence 
takes place and in the rates of convergence. Both questions strongly depend on m. Let us 
emphasize for instance that the Barenblatt solution Ud,t is integrable in y for m > m c , while 
the pseudo-Barenblatt solution corresponding to m < m c is not integrable. Since much is known 
in the case m > m c , see for instance [HJ [25] and pUl Ell EHJ EJ E51 E7J EU US] for more complete 
results, the main novelty of our paper is concerned with the lower range m < m c , which has 
several interesting new features. For instance, in the analysis in high space dimensions, that is 
d > 4, another critical exponent appears, 

d-4 
m*:= — <m c . 

A key property of is that the difference of two pseudo-Barenblatt solutions is integrable for 
m G (m^,m c ), while it is not integrable for m G (0,771*]. 

The convergence towards Barenblatt and pseudo-Barenblatt solutions is subtle since the solu- 
tions converge to zero everywhere. To capture the asymptotic profiles, it is therefore convenient 
to rescale the solutions and replace the study of intermediate asymptotics by the study of the 
convergence to stationary solutions in rescaled variables, 

d.3) t: = log (|sy) - d * = 

with R as above. In these new variables, if u is a solution to (jl.ip . the function 

v(t,x) :=R(T) d u(T,y) 
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solves a nonlinear Fokker-Planck type equation, 

( d t v(t, x) = Av m (t, x) + V ■ (x v(t, x)) (t, x) G (0, +00) x R d , 

(1-4) 

[ v(0,x) = v (x) x eR d . 

The initial data for (jl.lj) and for the rescaled equation (jl.4l) are related by 

Mv) =R(0)- d v o (y/R(0)) , 

where R(Q) = [d | m — m c \ T \ d(m_mc) only depends on T. In this formulation, the Barenblatt and 
pseudo-Barenblatt solutions are transformed into stationary solutions given by 

(1.5) Vd {x):=(d+ 1 —^\xA " m 

\ 2m J 

where < m < 1 and D > is a free parameter. With a straightforward abuse of language, we 
say that Vr> is a Barenblatt profile, including the case m < m c . The value D = can also be 
admitted as a limit case, but the corresponding solution is singular at x = 0. See [37] for more 
details. The parameter T has disappeared from the new problem, but it enters in the change of 
variables. Note that in all cases, t runs from to infinity in these rescaled variables. 

Assumptions and main results. We can write the assumptions on the initial conditions in 
terms of either uq or vq. We assume that 

(HI) uq is a non-negative function in Lli oc (M rf ) and that there exist positive constants T and 
Dq > D\ such that 

U Do , T (0, y) < u (y) < U Dl , T (0, y) V y G R d . 
(H2) There exist G [£>i, Ar] and / G Ll(R d ) such that 

Mv) = u D ., T fry) + f(y) VyG]R d . 

Note that by the Comparison Principle, see Lemma 12.21 b elow . in the case m < m c , (HI) implies 
that the extinction occurs at time T. When m > m*, (H2) follows from (HI) since the difference 
of two Barenblatt solutions is always integrable. For m < m*, (H2) is an additional restriction. 
We shall assume throughout this paper that d > 3 and observe that (H2) has to be taken into 
account only if m* > 0, that is, d > 5. 

In terms of vq, with / replaced by R(0)~ d f(y/R(0)), conditions (HI) and (H2) can be rewritten 
as follows. To avoid more notations, we keep using / in (H2') although it is not the same function 
as in (H2). 

(HI') vq is a non-negative function in Lli oc (R rf ) and there exist positive constants Do > D\ such 
that 

V Do (x)<v (x)<V Dl (x) V x G R d . 
(H2') There exist D* G [D^Dq] and / G Ll(R d ) such that 

v (x) = V Dt (x) + f(x) VxeR d . 

If m G (m*, 1), the map D 1— » J R d(^o — Vd) dx is continuous, monotone increasing. Hence we 
can also define a unique G [D± , Dq] such that 



/ («o - Vb.) dx = 



Before stating any result, one more exponent is needed. We define p(m) as the infimum of all 
positive real numbers p for which two Barenblatt profiles Vd 1 and Vb 2 are such that \Vd 1 — Vd 2 \ 
belongs to L p (R d ): 

P[ ' ' 2 (2 - m) ' 
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We see that p(m) > 1 if m G (0, m*), p(m*) = 1, and p(m) < 1 if m > m*. 

We can now state the convergence of v(t) towards a unique Barenblatt profile. For simplicity, 
we will write v(t) instead of x *— > v(t, x) whenever we want to emphasize the dependence in t. 

Theorem 1.1 (Convergence to the asymptotic profile). Let d > 3, m G (0,1). Consider the 
solution v of (|1.4p with initial data satisfying (H1')-(H2'). 

(i) For any m > m* ; there exists a unique G [Di,Dq] such that f^d(v(t) — Vd„) dx = 
for any t > 0. Moreover, for any p G {p(m), oo], lim^oo f Rd \v(t) — VdJ p dx = 0. 

(ii) For m < m*, v(t) — Vn* is integrable, f R d(v(t) — Vb*) dx = J" Rd / dx and v(t) converges 
to Vd* in L p (]R rf ) as t — > oo, for any p G (1, oo]. 

(iii) (Convergence in Relative Error) For any p G (d/2,oo], 

lim \\v{t)/V D .-l\\ p = 0. 

In case m > m*, the value of can be computed at t = as a consequence of the mass balance 
law f K d(vo — VdJ) dx = 0, and then the conservation result holds for all t > 0, see Proposition 12.31 
below. On the other hand, in the case m < the mass balance does not make sense, but -D* 
is determined by Assumption (H2). In this case, the presence of a perturbation / of Vd* with 
nonzero mass, does not affect the asymptotic behavior of the solution at first order. 

In a recent paper [22], P. Daskalopoulos and N. Sesum prove some of the results of Theorem ll.il 
under similar hypotheses (see [221 Theorem 1.4]). Actually they only prove the L°° convergence 
in case (ii) and the LI n L°° convergence in case (i). Our proof was obtained independently and 
announced in [7]. It is based on entropy estimates and paves the way to the sharper results on 
convergence with rates, which are the main purpose of the present paper. Assertion (iii) says 
that the convergence of (i)-(ii) can be improved into a convergence in relative error, in the sense 
of [15]. Such a strong convergence may look surprising at first sight, but it is a consequence of 
Assumption (HI'): the tails of vq and Vo t have the same behavior as \x\ — > oo. 

We can now state our main asymptotic result, on rates of convergence. To state this second 
result, we need yet another exponent, 

2d(l-m) 
q * '~ 2(2 - m) +d(l - m) 
and note that > 1 if m < m*, q* = 1 if m = m*, and q* < 1 if m > m*. For any q > g*, 
the function Vd* is in 

L (2-m) 9 /(2- 9 )( R< ^ which allows us tQ 

use convenient Holder interpolation 

inequalities. We define the C 3 semi- norm by 

II/IIcj(R*) := ™ ax SU P \ dT, f( x ) | > 

\m=j xem. d 

where the standard multi-index notation is used: |r/| = rji + . . . + r)d is the length of the multi- 
index i] = (771, ... , rjd) G The last ingredient is a Hardy-Poincare constant, which is defined 
as follows. For any m G (0,m») U (m*, 1), let 

(1.6) , mj:=roiDf k<WVn.4* 

>• /„ a |ft-fci 2 ^: m dx' 

where the infimum is taken over the set of smooth functions h such that supp(/i) C M. d \ {0} 
and h = if m < m*, while for m > m*, h := J Rd hV^~ m dx/ J Rd V^~ m dx, cf. Theorem IA.1I in 
the Appendix for more details. We shall prove tlictt A m d is positive and independent of D*. In 
the next result, the time decay rate is formulated in terms of the spectral gap X m d- Analyzing 
the relationship between the optimal constant C m ^d = m /^md m the corresponding functional 
inequality and the asymptotic rates of the fast diffusion equation is the leitmotiv of this paper. 
The case m = m* has to be excluded for reasons which are deeply related to Hardy's inequality, 
see [7]. 
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Theorem 1.2 (Convergence with rate). Under the assumptions of Theorem \l.l\ with \ m ^d given 
by (|1.6j) . if to m*, there exists to > snc/i £/ia£ £/ie following properties hold: 

(i) For are?/ q G , oo] ; i/iere exists a positive constant C q such that 

K<) -vb.|| g <c ?e - A ^* vt>t . 

(ii) For are?/ 1? G [0, (2 — m)/(l — m)), i/tere exists a positive constant K$ such that 

|| |a:|*(u(t) - Vb.)|| 2 < K 6 e~ x ^ dt Vt > t . 

(iii) For any j G N, i/iere exists a positive constant Hj such that 

2 A m,d . 

II « (*) - Vb. lb(R<*) < Hj e d + 2 ^'+ 1 ) V t > t • 



The constants C q , K$ and Hj depend on to, m, d, vq, Dq, D±, and q, $ and j; to also depends 
on Dq and D\. It is remarkable that the decay rate of the nonlinear problem is given exactly by 
A m ,d (see Section . Using (jl.3p . the results of Theorem 11.21 for the solution v(t) of (|1.4j) can 
be translated into results for the solution n(r) of (jl.lj) as follows. 

Corollary 1.3 (Intermediate asymptotics) . Let d > 3, m £ (0,1), m ^ m*. Consider a solu- 
tion u of (jl.ip with initial data satisfying (H1)-(H2). For r large enough, for any q G (o*,oo], 
there exists a positive constant C such that 

\\u{T)-U Dt {T)\\ q <CR{rr a , 

where a = X m! d + d(q — l)/q with \ m ,d given by (|1.6j) . and large means T — r > 0, small, if 
m < m c , and r — > oo if m > m c . 

We also obtain a convergence result in relative error. For any p G (<i/2, oo], define 

A(P) := f(d + 2)(2-m) Xm ' d ■ 

Theorem 1.4 (Exponential Decay of Relative Error). Under the assumptions of Theorem 11.21 
if m ^ m* ; /or any p G (d/2, oo], i/iere exists a positive constant C and A G (0, A(p)) suc/i that 

\\v(t)/V D , -l\\ p <Ce~ xt Vt>0. 

Let us list a few observations on the above results. 

(a) In the range m c < to < 1, convergence with rates has been obtained under various as- 
sumptions, cf. [25l HH [15] (optimal rates) for to G [toi,1), mi = (d—l)/d, and [TH [16] for 
to, G (?n c ,TOi), which are weaker than the ones of Theorem II. 2[ See [27] for the detailed analysis 
of the spectrum of the linearized operator in the range to > m c . A stronger convergence has also 
been proved in the sense of relative error under very mild assumptions, cf. |45j . 

(b) In the range m c < to < 1, Assumption (HI) is less restrictive than one could think. By the 
global Harnack principle of [11], see Theorem 12.51 below, any solution with non-negative initial 
data no G Lli oc (M rf ) that decays at infinity like uo(y) = Odyl 2 ^ 1 ""^), is indeed trapped for all 
t > between two Barenblatt solutions if m c < to < 1. The restrictions on the class of initial data 
are therefore not so essential as far as the asymptotic behavior is concerned, and can therefore 
be relaxed. 

(c) In the range < to < m c , the pseudo-Barenblatt solutions are not integrable. For m < m c 
many solutions vanish in finite time and have various asymptotic behaviors depending on the 
initial data. Solutions with bounded and integrable initial data are described by self-similar 
solutions with so-called anomalous exponents, see [341141] and [471 Chapter 7). Even for solutions 
with initial data not so far from a pseudo-Barenblatt solution, the asymptotic behavior may 
significantly differ from the behavior of a pseudo-Barenblatt solution: in |22[ Theorem 1.4] a 
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solution of (jl.lh is found, which is such that lim^^^ uq(x)/Ud,t(0, x) = 1 and which, after 
rescaling, does not converge to Vd as t — > oo, that is for r — > T. Assumption (H1)-(H2) are 
therefore more restrictive than for m > m c . 

(d) Proofs are constructive and the values of the various constants are explicit although not so 
easy to write. The interested reader will be able to recover their expressions by carefully reading 
the proofs, where all details are given. See Appendix A for more details on the constant X m .d that 
controls the rate of convergence. The rate given by this exponent is sharp in the linear case, and 
a deeper analysis should prove that it is sharp also in the nonlinear case. Obtaining the optimal 
value of \m,d is still an open question in the range m G (m^,m c ]. Our method gives convergence 
with rates even in the limit case m = m c , which is new. 

Further comments. After stating our main results, let us come back a little bit on the mo- 
tivations of this paper, on the main tools and the originality of our results with respect to the 
existing literature. 

During the last few years, asymptotic rates of convergence for the solutions of nonlinear diffu- 
sion equations have attracted lots of attention, usually in connection with time-dependent scalings 
and entropy methods. This has been first done in the range of exponents corresponding to the 
porous medium equation, with 1 < m < 2, and in the range where standard Gagliardo-Nirenberg 
inequalities apply, mi < m < 1, see [251 El CE] • The class of non- negative, finite mass solutions 
has to be narrowed to the smaller set of functions with finite free energy, or to be precise, with 
finite entropy and finite potential energy. In the rescaled variables, asymptotic stabilization to 
the Barenblatt profiles holds at an exponential rate, while in the original time variable r, the 
convergence of the difference with the Barenblatt solutions holds at a power-law rate, which is 
shown to be optimal. 

The next question was to understand what happens for m < mi. After the linearized analysis 
of [EL the proof of convergence with rates was done in [16] in the range m c < m < mi for which 
global existence of finite mass solutions still holds. The basin of attraction is narrowed to the 
class of solutions with finite relative entropy with respect to some Barenblatt solution. 

A dramatic change occurs for m < m c , since a large class of solutions vanish in finite time. 
As a consequence, mass is not conserved, and a key estimate for higher values of m is lost. It 
is however natural to investigate the basin of attraction of the pseudo-Barenblatt solutions for 
m < m c using relative entropy techniques and to study the convergence rates. This can be done 
in a weighted space using functional inequalities, which can still be related to some spectral 
properties of a differential operator obtained by an appropriate linearization. 

The generalized entropy functional, or free energy functional, is defined as 



It is then observed that the free energy of the Barenblatt profiles, cf. \15\ 125]. becomes infinite if 
m < m o, where := d/(d + 2) G (m c ,mi). In order to avoid this difficulty, it is convenient to 
work with the relative entropy of v with respect to Vd defined as follows: 



The relative entropy is the key tool of our analysis. It is such that £[u|Vd] := £[v) — £[Vd] if 
m G (mo, 1) and L d v dx = J Rd Vd dx, that is for D = D*, where is as in Theorem I l.li The 
functional £[i>|VdJ can also be defined for m < mo. By homogeneity of (p, we can indeed rewrite 
it as 




where 




£[v\V D ] : 




ip(v) - tp{V D ) - ^'(Vb) (v - V D ) dx . 




<p(w) - ip(l) - tp'(l) (w - 1) VSl dx with w 



'V 
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This makes clear why it is well defined at least for w close enough to 1 as |x| — ► oo. The functional 
v i — ► ^[v|yb,] is convex and achieves its minimum, 0, for v = Vd»- If u is a solution of (|1.4j) . the 
entropy production term takes the form 

-±£[ v (t)\V D ,]=l[v(t)\V D ,], 



where the functional 

v ^ l[v\V D ) := I v V(f'(v) - V(f/(V D ) 



2 

dx 

will be called the relative Fisher information. See Proposition 12.61 for more details. For any 
m G [mi, 1), £MVb„] < ^Z[u|VdJ holds for any smooth function v and the inequality is nothing 
else than the optimal Gagliardo-Nirenberg inequality, for which equality is achieved precisely by 
the Barenblatt profiles, see [25]. In such a case, 

S[v(t)\V D ,]<£[v \V D ,]e~ 2t Vt>0. 

The limit case m = mi corresponds to the critical Sobolev inequality whose optimal form was 
established by T. Aubin and G. Talenti in [HE], while in the limit m — > 1 one recovers Gross' 
logarithmic Sobolev inequality, see [29U25] . For m G [mi, 1), F. Otto in [3D] noticed that (jl.4l) can 
be interpreted as the gradient flow of the free energy with respect to the Wasserstein distance. 
The exponent m = m\ is the limit case for which the displacement convexity property holds true. 

Pushing the method to the case < m < mi requires the use of the relative entropy in place 
of the free energy. The method applies only to a class of initial data which have a finite relative 
entropy with respect to some Barenblatt profile Vd* and satisfy convenient bounds. Mass can 
be finite in the case m G (m c , mi), which was the framework of some earlier studies, see |14[ I16j. 
or infinite if m G (0, m c ). Two Barenblatt profiles Vd and Vd 1 have finite relative entropy, i.e. 
£[Kdi|Vd ] < oo if and only if either d < 4, or d > 5 and m > m*, = (d — 4) / (<i — 2). Hence, 
for d > 5, m = m* is a threshold not only for defining the relative mass of two pseudo-Barenblatt 
solutions, but also for defining their relative entropies or for the integrability of V^~ m . Note 
that m* < m c for all d > 5. The proof of Theorem 11.21 amounts to prove that the relative 
entropy £[v|VdJ decays in time and converges to at an exponential rate when t — > oo. For 
m > min{0, m*}, £[u|Vb„] is well defined under condition (HI'). For m < m*, an additional 
restriction is required, which is precisely the purpose of (H2'). 

Our approach of course covers the case m > m c and we recover some of the results found 
in |14[ll6j. Some of our results can also be extended to the range m < 0, but additional technical 
complications arise, which are still to be studied. In this paper, we leave apart several interesting 
questions, like the precise study of the case of m = m* or the equation ut = Alogu in dimension 
d > 2, see e.g. |20[ [21] [42l 148] . which is the natural limiting equation to study in the limit 
m — ► 0. Also see [HJ [32], [33] for results which seem closely related to ours, and [26] in the 
case m = (d — 2)/(d + 2). In particular we do not use the Bakry-Emery method introduced in 
[2], on which the results of \15\ [13] [36l [14] [16] are based. We prove a conservation of relative 
mass, which allows us to remove the limitation m > m c . Neither mass transportation techniques 
nor Wasserstein distance are needed, although the approach of Section [A. 31 is not unrelated, see 

EIEIHEH]. 

The paper is organized as follows. In Section [21 we extend the property of mass conservation, 
which holds only for m > m c , to a property of conservation of relative mass, see Proposition 12.31 
This selects a unique Barenblatt profile, which governs the asymptotic behavior. We also establish 
regularity properties of the solutions. From there on we work with the quotient of the solution 
by the Barenblatt profile. In Section we prove Theorem 11.11 and establish several properties 
which are used in the sequel. Sections [4] and [5] are respectively devoted to introducting a suitable 
linearization and to the derivation of entropy - entropy production estimates in the nonlinear 
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case, from the corresponding spectral gap estimates for the linearized problem. The proof of 
Theorem 11.21 is given in Sectional 

Appendix A is devoted to the proofs of spectral gap estimates, that is, of weighted Poincare- 
Hardy inequalities, which have already been studied in [7J, and in [27] in the special case 
m > m c . We consider the family of weights of the form Vrj or V^~ m , D > 0, that are ob- 
tained from the linearization of the relative entropy. In the limit D — > 0, they yield the case 
corresponding to the weighted L2 norm of the Caffarelli-Kohn-Nirenberg inequalities, cf. [121 117] . 

A final section, Appendix B, explains how to extend the results of this paper to the fast 
diffusion with exponents m < 0. Note that the equation needs to be properly modified. The 
conclusion is that the results still hold and the proofs need only minor modifications that are 
indicated. 

2. Basic Estimates 

We establish in this section the main properties of the solutions that will be used in the sequel. 

2.1. Ll-contraction and Maximum Principle. 

Lemma 2.1 (Ll-contraction). For any two non-negative solutions u\ and u 2 of (11.11) defined 
on a time interval [0, T), with initial data in Lli oc (IR rf ) ; and any two times t\ and t 2 such that 
< t\ < t 2 < T, we have 

\ \u\{t 2 ) - u 2 (t 2 )\ dx < \ \u x (ti) - u 2 (*i)| dx . 

JR d JR d 

The above result is well-known to be true for solutions with LI data, cf. Proposition 9.1], 
even in the stronger form 

/ \ui(t 2 ) - u 2 {t 2 )\ , dx < / [ui(h) -u 2 {h)] dx , 

JR d JR d 

where [u]+ denotes the positive part of u. The result for data in Lli oc (IR d ) is obtained by approx- 
imation, using the uniqueness of solutions to the Cauchy Problem, which has been established 
in [30]. Note that when the right-hand side is infinite the result applies but there is nothing to 
prove. As a consequence, we also obtain the following. 

Lemma 2.2 (Comparison Principle). For any two non-negative solutions U\ and u 2 of (jl.ip 

on [0, T), T > 0, with initial data satisfying uq\ < uq 2 a.e, uq 2 G Lli oc (]R rf ) ; then we have 
ui(t) < u 2 (t) for almost every t G [0, T). 

We will see below that under Assumption (H1)-(H2) the solutions are smooth functions, hence 
the comparison in the previous result holds everywhere in [0, T) x W 1 . 

2.2. Conservation of relative mass. Mass conservation is used in the range m > m c to 
determine the parameter D which characterizes the Barenblatt profile Vd- In the range m < m c , 
we can still prove that J R d(v(t) — Vd) dx is conserved for any t > 0, even if Vb Ll(M rf ). 

Proposition 2.3. Let m G (0, 1). Consider a solution u of (11.11) with initial data uq satisfying 
(H1)-(H2). If for some D > 0, J K d(wo — Ud,t(®, •)) dx is finite, then 

[ [u(t, x) - U DiT (t, x)] dx = f [u (x) - U DiT (0, x)] dx V r G (0, T) . 

jR d JR d 

Proof. In the range m > m c , uq and Ud,t(0, are integrable and mass conservation of the 
solutions of (jl.ip is well-known. 
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Assume next that to < m c and let x be a C2 function on R + such that % = 1 on [0, 1], x — 
on [2, oo), and < x < 1 on [1) 2]. For any A > 0, take 4>\{x) := x(\x\/X) as a test function and 
denote by B\ the ball B(0, A). Then, 



I [u(r) - U d , t (t)} <j> x dy = f [u m (r) - U^ T {r)} A0 A dy 
dr J R d J R d L 



B 2 \\B\ 

< c 



[n m (r) - U^ t {t)] A0 a dy 

\u(t)-U d>t (t)\ USJ(t) |A0 a | dy, 

B 2 \\Bx 

where C is a numerical constant depending on D, Dq, D\. As A — > oo, we observe that in 
-E>2A \ -B A , Up~ T and |A0 A | behave like A2 and A~ 2 , so that Uq~ t (t) A0 is bounded uniformly 
with respect to A. The right hand side is therefore bounded by f B2X \ Bx \ U ( T ) ~ Ud,t(t)\ dy. For 
any n, T2 G [0, T), we write 



KT2) - ?7d,t(t2)] 0a dy 



[u(n) - Ud s t(ti)] 4>\ dy 



T2 r 

< I / \u — Ud,t\ dy dr 
'n Jb 2X \b x 



By the Ll-contr action principle, see Lemma 12.11 we also know that |u(r, y) — Uo,T(j,y)\ is 
integrable in y, uniformly for all r > 0. The integrability condition implies that the right-hand 
side goes to zero in the limit A — ► +00. 

The case m = m c is similar, except that there is no extinction time. □ 

In the rescaled variables given by (jl.3p . relative mass is also conserved. Consider a solution v 
of (|1.4p with initial data vq satisfying (H1')-(H2'). If for some D > 0, J R d(vo(x) — Vb(i, %)) dx is 
finite, then 

/ [v(t, x) - V D (t, x)] dx= [v (x) - V D {t, x)] dx V t > . 

Whenever m < m c , we recall that pseudo-Barenblatt solutions have infinite mass, that is 
J Rd Vd„ dy = 00, but we observe that the difference of two pseudo-Barenblatt solutions is inte- 
grable if m > m*. In such a case, the parameter D G [Di, Do] in Proposition 12.31 can be arbitrary. 
In the proof, we can moreover estimate f B2X \ Bx W — U b ,t\ dy by 



max ■ 



B 2\/R(t)\B\/R(t) 



\V Do -V D \dy, f \V Dl -V D \dy 

J B 2 \/R(t)\Bx/R(t) 



which converges to as A — > 00. As already quoted in the introduction, the map D 1— > f Rd (vo — 
Vd) dx is continuous, monotone increasing, and we can define a unique G [D\, Dq] such that 
JrodC^O — Vb,) dx = 0. Under the assumptions of Proposition 12.31 



[v(t,x)-V Dt (x)]dx = v*>o. 

This fact is used in the statement and proof of Theorem 11.11 for to > to* . On the contrary, if 
m G (0,771*], integrals are infinite unless D = D* in Proposition 12.31 and then, with the notations 
of Assumption (H2'), 



[v(t,x) - V Dt (x)] dx 



f dx V t > . 



In both cases, that is for any to G (0, 1), we shall summarize the fact that — Vb„] dx = 

by saying that the relative mass is conserved. 
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2.3. Passing to the quotient. Consider a solution v of (|1.4h . As in (45j CEH [H3] , define 
(2.1) 



Vb.(x) 

Next, we rewrite Problem (jl.4p in terms of w. 
1 



w (t, x) := V (t, x) E (0, oo) x R d . 



(2.2) 



Define 



^vf^- f («; m - 1 -l)VS:- lN ) 
\m — 1 * / 



in (0, +oo) x 



w(0, •) = -u; - 



_vo_ 

Vd, 



Wq := inf 



m 



V 



Do 



< 



sup 



Vd, 
Vd, 



:= Wi . 



A straightforward calculation gives 

Wn = — 



In terms of i«o, assumptions (HI') and (H2') can be rewritten as follows: 

(HI") wq is a non-negative function in Lli oc (R rf ) and there exist positive constants -Do > D\ such 
that 

< W < ^44 < w{x) < ^44 < W 1 < +oo V x E R d . 



(H2") There exists / E Ll(M d ) such that 

w(x) = 1 + 



VdM 

m 



Vie 



VdM 

As a consequence of the Comparison Principle, see Lemma (221 (HI") is satisfied by a solution w 
of (|2.2p if it is satisfied by wq- 



2.4. Regularity estimates and Harnack principle. We start by briefly recalling some well- 
known results for solutions of fast diffusion equations. A basic regularity result is due to 
DiBenedetto et al. , see |18|. p. 270], and concerns local space-time Holder regularity for Prob- 
lem (jl.ip . with some Holder exponent a E (0, 1); it holds for locally bounded initial data, possibly 
with sign changes. In the present situation of locally bounded and positive initial data, it is known 
that the solutions are C°° as long as they do not vanish identically because we avoid any degen- 
eracy or singularity and the standard parabolic theory applies. In the sequel we are in particular 
interested in some sort of uniform CI regularity under the assumption (HI"). We find that it 
is preferable to work with the function w introduced in (|2.ip since it is uniformly bounded from 
above, and from below away from zero.- 

Theorem 2.4 (Uniform C k regularity). Let m E (0, 1) and w E L^ c ((0,T) x R d ) be a solution 
of ([22]). Then for any k E N, for any t E (0, T), 

sup \\w(t)\\ C k iRd) < +oo . 

t>t 

Proof. Take t > to > 0. For a given A > 0, the equation for w is uniformly parabolic in B\, so 
we conclude that the regularity estimates hold on for any t > to- Let v be the solution of 
(ll.4p corresponding to w. In order to cover the large values of we consider the scaling 

v x (x,t) = A 2/(1 - m) t>(Ax,t) 
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with A — > oo. Then v\ is again a solution of (|1.4j) . but the region ft A = {x G R d : A < |z| < 2 A} 
gets mapped into the region $7i = {x G R d : 1 < |x| < 2}, for all A. Note also that this scaling 
transforms the Barenblatt profiles according to 

(Vd)a = V D /\ 2 , 

so that, on Qi, (Vd)a is uniformly bounded from above and from below in Qi as A — > oo, and 
converges to Vo- 

Next, we pass to the functions w\(x,t) = v\(x,t)/(VD,)\(x) and observe that in fii, (wa)a>i 
is uniformly bounded from above, and from below away from zero. Since w\ satisfies (|2.2p with 
Vd* replaced by (Vd»)a> we conclude as in part (i) that 

\\wx{t)\\ ck{p/l) < C k 

uniformly in t > to and A > 1. Undoing the scaling, we find a constant C independent of A such 
that, for any t > to, 

\Vw{t,Xx)\ C 
w(t,Xx) ~ A 

We conclude that the result holds for k = 1. The same argument applies for k > 1. □ 

In the range m G (m c , 1), the regularization effects of the fast diffusion equation allow to get 
rid of assumption (HI). See [11] for details expressions of the constants entering the statement. 

Theorem 2.5 (Global Harnack principle, Theorem 1.2]). Let uq be a non-negative function 
in Ll(R d ) solution of (jl.ip . Assume that for some R > 0, sup| y | >i? uo(y) \y\ 2 /( l ~ m ) is finite. 
Then, for any e > 0, there exist positive constants To, T\, Do and D\, such that 

Ud ,t (t, y) < u(r, y) < U Di>Ti (t, y) V (r, y) G (e, oo) x R d . 
2.5. Relative entropy. In terms of w, we define the relative entropy 

(2.3) m := 7—— / [(w-l)-±-(w m - 1)] VS dx 

i — m j^d m 

and the relative Fisher information 

i2A) j w : = (^1)2 1 1 v - vsr'i l 2 » d, . 

These definitions are consistent with the ones given in the introduction, in the sense that, for 
w = v/ Vd, , 

F[ w ] = —£[v\V D ] and J[w] = —X[v\V D ] . 
m m 

The 1/m factor simplifies the expressions of the linearized relative entropy and Fisher information, 
as we shall see below. It has no impact on the rates. Consistently with the passage to the quotient, 
the relative entropy and the relative Fisher information are related as follows. 

Proposition 2.6. Under Assumptions (H1")-(H2"), if w is a solution of (12. 2j) . then 
(2.5) ^ nw{t)] = - J[w{t )) . 

Proof. As in Section \2. 2\ consider a test function 4>\(x) := x(|x|/A) where \ is a smooth function 
on R+ such that x = 1 on [0, 1], x = on [2, 00), and < x < 1 on [1, 2]. Then, using ([22]), the 
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equality (w 



m—l 



DV, 



m—l „,m— 1 



D, 



V 



1 and integration by parts we obtain 



— !— / [{w-l)--{w m -l)]4>xV^dx 
1 - m J R d m 



d 1 



m 


(m 


-1)2 




m 


(m 


-1)2 




m 


(m 


-1)2 



V^™" 1 _ i) vg*-i] . vf^" 1 " 1 - 1) V™" 1 A ] wV D , dx 



f V (t- m - 1 - v™~ 1 ) • [^v^- 1 - y^- 1 ) + vMv™- 1 - v™- 1 )] v dx 



where the last integral can be computed as 



1Z{\) :-- 



m 



2{m- 1)2 
m 

' 2(m - 1)2 



,m— 1 T/m— 1 



,m— 1 T/m— 1 



I?, 



(Vv • V(/)\ + w • A0 A ) dx 



In the region S7a = {x G R d : A < |x| < 2 A}, A > 0, we get 

\K{\)\ < h I \v- V D f V^~ 2) (\Vv\\V<f> x \+v\A<f> x \) dx 



< ki sup 



v 2(rn-2) ^ _ y ^ ( | Vu | | V ^ A | + t> | A^ A | ) 



\v — Vb* | dx , 



where the positive constant k\ depends on m, d, Dq and D\. 

As in the proof of Theorem 12.41 consider a solution v of (jl.4p and define 

v x {x,t) = \ 2/{1 - m) v(Xx,t) . 

In what follows, c% will denote positive constants which may depend on m, d, Dq, D\ and on the 
maximum of V</>i, but not on A. 

For any A > 0, v\ is a solution of (|1.4p but the region Q\ = {x G R : A < |x| < 2 A} gets 
mapped into the region = {x S M d : 1 < |x| < 2}. We already know that Vv\ is uniformly 
bounded on fii, by Theorem 12. 4t sup^g^j V«a(*, #)| < 0)- In terms of v, this gives the estimate 

A i-™ sup |V y f(t,y)| = A 1 -™ sup |A _1 V a .'y(t, Ax)| = A -1 sup |V x UA(t, x )| < A -1 c 

y££l\ x£Qi x£Qi 

and proves that 

By our choice of 4>\ , we see that 



sup |Vt>(£,x)| < Co A 



■ -i 



sup|V0 A |<-r and sup|A</> A |< — 



Putting together these two estimates, we get 



|Vu| |V0a| +v\A(f> x \ < c 3 A i- 



Next, we observe that 
dV D 



dD 1 — m 



D + 



1 — m 
2 m 



2-m 
1 — m 



1 



1 — m 



2-m 
D 



Hence, for some constant C4 depending on m, d, Dq and D\, 



\V Dl - V Do \ <c 4 ^ 



2-m 
Do 
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and 



sup 



2(m-2) 



Dn 



\v- V D .\ (\Vv\\V<p x \+v\A(b x \) 



< C4 sup 



V^- 2 (|Vt>||V^ a |+t;|A^ 



Taking into account the fact that, for any A > 0, Vd < C5 A 2 ^ 1 m ) on we obtain 



sup 



2(m-2) 



D 



v-V Dt \ (\Vv\ \V(j> x \ + v|A0 A | ) 



o 2 — m 2 o 

< C6 A i-" 1 A !- m = C6 



for some positive constant cq which is independent of A. By assumptions (H1')-(H2') and the 
Ll-contraction principle, the difference v — Vd, is in LI, and so, lim^oo Jq \v — Vd,\ dx = 0. 
This proves that lim^oo 7£(A) = and we conclude by passing to the limit as A — > 00. □ 



3. Convergence without rate and in relative error 
This section is mostly devoted to the proof of Theorem 11.11 

3.1. Relative entropy. Under Assumptions (H1")-(H2"), the relative entropy is well defined. 
Lemma 3.1 (An equivalence result). Let m E (0,1). If w satisfies (H1")-(H2"), then 

\w™~ 2 [ \w-l\ 2 V% dx<F[w] < \w™~ 2 [ \w-l\ 2 VEldx. 

2 jRd 2 J R d 

Proof. For a > 0, let 4> a (w) := j^— [(w — 1) — (w m — l)/m] — a(w — 1) 2. We compute 4>' a (w) = 
[l - w™- 1 ] - 2a (w - 1) and <f/£(w) = w m ~ 2 -2 a, and note that <p a {l) = <j>' a (l) = 0. With 
a = W™~ 2 /2, (f>" is positive on (Wo, Wi), which proves the lower bound after multiplying by 
and integrating over R d . With a = W™~ 2 /2, <\>" a is negative on (Wo, Wi) which proves the upper 
bound. □ 

Lemma 3.2 (Boundedness of the free energy). Let m G (0, 1). If wq satisfies (H1")-(H2"), then 
the free energy J-[w(t)] is finite for any t > 0. 

Proof. By virtue of Proposition 12.61 we have to prove the result only for w = wq . Notice that 

q/ 2 — m \ 

for any Dq, D\ > there exist a positive constant C such that \Vd — VbJ < C \x\ y t-™> as 
\x\ — ► co. Indeed, 



\Vd ~V Di \ 



1 — m\ 1 ~ m 2m\Do — D 



^| x |- 2 (f^7)(l +0 (l)) 



as x — > 00 . 



2m J (l-m)2 

By Lemma l3.1( for some positive constant c depending on Dq and D\, we have 

2 

< [ ^- - 1 VR dx < [ \v - V D f V^ 2 dx < c [ \Vdo- V Dl \ 2 V™- 2 dx . 
Jm.d Vd, Jrci j R d 



in-: 



If m G (m*, 1), then |Vd — VbJ ^iT" = 0[\x\ i-mj i s integrable as |x| — ► 00. Otherwise, if 
m G (0,771*], F[w\ is also integrable as \x\ — ► 00 because 



W, 



Li < Tfi -?>o] < / |/| 2 V™- 2 dx < f I/) IVdo - V Dl I V^" 2 drr 







/ is integrable and |Vd — Vd 1 | V^ 1 is bounded (we ask the reader to check this fact). □ 
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3.2. Pointwise convergence in relative error. 

Lemma 3.3. Let m G (0,1). If w is a solution of (|2.2p with initial data wq satisfying (Hl")- 
(H2") ; then lim^oo w(t, x) = 1 for any x <E~R d . 

Proof. Let w T (t,x) = w(t + r, x). By the uniform C k regularity, see Theorem 12.41 the functions 
w T are uniformly CI continuous. Hence, by the Ascoli-Arzela theorem, there exists a sequence 
T n — > oo such that w Tn converges to a function w^, locally uniformly in (t, x). We know by the 
Comparison Principle, see Lemma f2.2t that Woo > Wq > 0. By interior regularity of the solutions, 
the derivatives also converge everywhere. 

By Lemma 13.21 T[w\ is finite. Since 

F[w(T n )] - T[w{r n + 1)] = / J[w(s)] ds = / J[w{s + r n )} ds , 



as a function of t, iJ[w Tn ] is integrable on [0, 1] and converges to zero as n — > oo, 

lim f f |V [{w?-\t,x) - 1) V%-\x)]\ 2 w Tn (t,x)V D ,(x) dxdt = 0. 



By Fatou's lemma, Woo = lim^ooU;^ satisfies V [{w^~ l - l) V^" 1 ] = a.e. in (0, 1) x R d . As 
a consequence of the conservation of relative mass, see Proposition 12.31 w oo = L Thus, we have 
proved the convergence a.e., and by equi-continuity, the pointwise convergence. Since the limit 
is unique, the whole family {w T } T converges everywhere as r — > oo. □ 

3.3. Proof of Theorem 11.11 

Proof of Theorem \l.ll (i) and (ii). By Lemma 13.31 lim^oo \v(t,x) — Vd*(x)\ = for any x € M. d . 
Moreover, we observe that 



\v(t) - V D ,\ < max{ \V Do - V D ,\ , \V Dl - V D ,\ } = O (\x 



-2(2-m)/(l-m) 

as \x\ — > oo. By Lebesgue's dominated convergence theorem, v(t) converges to Vb, in L p (R d ), for 
any p G (p(m), oo), where pirn) := 2(2-m) ^ s ^ ne hifimum of all positive p such that the difference 
between two different Barenblatt profiles belongs to L p (M d ). 

The uniform convergence is based on the following interpolation lemma, due to Nirenberg, 
cf. [391 P- 126]. Let A, \x and v be such that — oo < A < \i < v < oo. Then there exists a positive 
constant C\ t ^ tV such that 

(3-1) ||/|ir/- A < C X ,,, V ||/||^ A V / G C(R d ) , 

where || • \\i/ a stands for the following quantities: 

(i) If a > 0, then ||/|| 1/(T = (j Rd \f\^ dxf . 

(ii) If a < 0, let k be the integer part of (—ad) and a = \a\ d — k be the fractional (positive) 
part of a. Using the standard multi-index notation, where \r)\ = r/i + ... + % is the length of 
the multi-index r/ = (771, . . . %) G Z d , we define 



l/a 



\&i f{x) - &> f(y)\ 
max cff = max sup ; ; =: f rr«crodi if a > 

m=* 1 la m=*x,»^ \ x ~y\ a c } 

max sup \d v f(z) \ := ||/||c*(K d ) if a = 

I \m=k 



As a special case, we observe that = 117 \\ci{ 

(iii) By convention, we note ||/||i/ = sup z6H d \f(z)\ = ||/|| C o(Rd) 
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Let j e N and A = -(j + l)/d < fi = -j/d < v = 1/2 so that k = j + 1 and a = 0. 
Inequality (|3.1j) becomes 

2d d+2j 2 

( 3 - 2 ) ll/llcJ(IR d ) < C -(j+l)/d,-j/d,X/2 ll/llc^+i(R d ) "^" 2 

for any j G N. By applying this interpolation inequality f = v(t) — Vd* with j = 0, we obtain 

2d d 2 

(3.3) \\v(t) - VbJoo < q di „, 1/2 IK*) - VdJI^ |Kt) - Vb.||^- a • 

By Theorem 12.41 the CI norm is uniformly bounded. If q{m) < 2, that is, if d < 8, or d > 9 
and m > (d — 8)/(d — 4), we already know that lim^co \\v(t) — VbJ|2 = 0. Otherwise, we can 
interpolate \\v(t) — Vd„ H2 between ||u(t) — VdJIi < ll^o - Vb«||i (see Lemma [27T]) and \\v(t) — Vb» ||<j 
for some </ > q(m). This proves that lim^oo ||u(t) — Vb» ||oo = 0. □ 



Proof of Theorem \l.l\ (Hi). 

Corollary 3.4 (Uniform convergence of the relative error). Let m G (0,1). If w is a solution 
of (12. 2j) tuii/i initial data wq satisfying (H1")-(H2") ; i/ien 

lim — l||oo = . 

t^oo 

Proof. Because of the convergence of v(t) to Vd, m L°°(]R' : '), we know that w(t) converges uni- 
formly to 1 on any compact set of M d . By Assumption (HI'), v(t) is sandwiched between two 
Barenblatt profiles that have the same asymptotic behavior when \x\ is large. In terms of w, 
this means that \w(t,x) — 1| is small for \x\ large, uniformly in t. Global uniform convergence 
follows. □ 

The fact that w(t) converges uniformly to 1 as t — > 00 allows us to improve the lower and 
upper bounds W^o and W\ for the function w(t), at the price of waiting some time. For any e > 
there exists a time to = *o( £ ) > such that 

1 - e < w(t, x) < 1 + e V (t, x) G (t , oo) x R d . 

Corollary 3.5 (L p Convergence of the relative error). Let m G (0, 1). If w is a solution of (|2.2p 
with initial data wq satisfying (H1")-(H2"), then w(t) converges to 1 in L p (M. d ) for any p G 
(d/2, 00]. 



Proof. By Assumptions (H1")-(H2"), there exists a positive constant cq such that wo — 1 is 
bounded and for \x\ large, 



\w - 1| 



V d , 



< Vbi - Vdq < c 



Vb. 



1 + \x\ 



By Lemma |2.2| the same estimate holds for w(t). Hence w(t) - 1 £ L' for any q > d/2. Let 
5= (p — d/2) /2 > 0. By Holder's inequality, 



lim / \w(t) — l\ p dx < lim \\w(t) — 1]]^ / ( — 
Theorem 1 1.1( (iii), results from Corollaries 13.41 and 13.51 



Co 



+ \ x \ 



S+d/2 



dx = . 



□ 
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3.4. Uniform convergence and C a regularity. 

Lemma 3.6. Let m G (0, 1). Consider a solution v of (11. 4j) with initial data vq satisfying (Hl')- 
(H2'). There exists to > 0, a G (0, 1) and a positive constant 7i such that h(t) := v(t) — Vd* is 
in C a and 

(3.4) \\h(t)\\c°Qi*) < ^IIM*)IU Vt >t . 

Proof. Since both v and Vd* are solutions to equation (jl.4p . /i solves 

^ = V • [m(Vb, + /i) m_1 V /i + m ((Vb. + ^) m_1 - V^ _1 - V™~ 2 /i) v Vb.] • 

Let A > 0. By Theorem ll.il we know that for some to > 0, for any t > to, ||/i(t)||oo can be taken 
uniformly small and v uniformly positive on i?2A- We apply the standard quasilinear parabolic 
theory, see e.g. [351 Theorem 1.1, p. 418], with structure functions a,i(x,t,h,£) = A£ + Bh 

and a = 0, where A(x,t) := my™' 1 , B(x,t) := m [( — v _ y °* Vg~ 2 )W D ,]. Hence there 

exists a Holder exponent a G (0, 1) and a constant depending on the uniform bounds for the 
coefficients, and on A, such that (|3.4[) is verified in B\ x (to + l,oo). To extend the estimate 
uniformly to the whole space, x G M d , we use the same scaling argument as in the proof of 
Theorem 12.41 We leave the details to the reader. □ 



4. Linearization 

In order to better understand the asymptotic behavior of the solutions of (|2.2p . we linearize 
the equation around the equilibrium, introducing a convenient weight. Let g be such that 



(4.1) w(t,x) = l + e ' Vt>0, VxGM d , 

for some e > 0, small. Substituting this expression in Equation (|2.2|) and letting e — > 0, we 
formally obtain a linear equation for g, 

(4.2) & = A m g where A m g := m V™~ 2 (x) V • [V Dm Vg] . 

The linear operator A m : L2(R d , V^~ m dx) -> L2(M d ,V^; m dx) is the positive self-adjoint 
operator associated to the closure of the quadratic form defined for <fi G C£°(M. d ) by 



(4.3) l[0] :=m / |Vg| 2 Vb* dx . 

See |23t Theorem 2.6] for more details. 

With the same heuristics, we linearize the relative entropy T and the relative Fisher infor- 
mation J , which provides the functionals F and I, where I is given by (|4.3p and F is defined 
by 

(4-4) F[g]:=l [ \g\ 2 V^ dx . 

* JR d 

Note that F[g] is the L2(R d , V%~ m dx)-norm up to a factor 1/2. If g is a solution of p~2l) . then 

(4-5) ^F[ 5 (t)] = -Ib(t)]. 

Still at a formal level, the conservation of relative mass amounts to require 
/ («o-Vb.)dx=/ (w-l)Vb. dx = e / tfV^dx 

JR d JR d JR d 
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in the limit e — > 0. Hence, it makes sense to require that J Rd g m dx = and use the spectral 
gap estimate, see (7j and Theorem IA.1I With C m ^ = raj\ m ^ we obtain 



(4.6) 



2F[g] < 



m 



which gives, for the solution of (|4.ip . an exponential decay of the relative entropy, 

F[s(t)] < e- 2A "- d * F[g(0)} Vt>0. 



In Sections [5] and [6l we will compare the relative entropy estimates for the solutions of (|2.2p 
with the ones of the linearized problem. This is the main ingredient of the proof of Theorem ll.2i 

The connection with the Fokker-Planck equation is easy to understand at the level of the 
linearized problem. In the limit m — > 1, we observe that 

U|2 



lim 

m—*l 



l/(l-m) 



Vb. = {2TrD*) d/2 n with fi(x) 



e 2D, 



Equation (|4.2p formally converges to the Ornstein-Uhlenbeck equation, 

gt = fi~ 1 V- (/iV 5 ) • 

The spectral gap inequality (|4.6p corresponds in such a limit to the well-known Poincare inequality 
with gaussian weight, 



d/i < / |V0|^ d^. V € C°°(R d ) such that 



dn = , 



where d// := /i dx. Note that in the Gaussian case, a logarithmic Sobolev inequality holds, 
see [29], 

.JU2 



/ M 2 iog(- 



do, < 2 



/ |V0| ; 



J** \<t>\ 2 dn, 

which is stronger than the Gaussian Poincare inequality. This is not the case with the mea- 
sure Vb» dx. Although the spectral gap inequality (14. 6p holds true, there is no corresponding 
logarithmic Sobolev inequality. 



5. More on the relative Fisher information 

In this section, we relate the relative Fisher and linearized Fisher informations. This and 
Lemma 13.11 provide us with an estimate of the relative entropy in terms of the relative Fisher 
information, or entropy - entropy production inequality, for the nonlinear problem. 

5.1. Fisher information and linearized Fisher information. 

Lemma 5.1 (Upper bound on the Fisher information). Let m G (0, 1). There exists two positive 
constants Pi and P 2 (depending on Wq, W\ and m) such that, for any w satisfying (H1")-(H2"), 

\[g]<PiJ[w]+P 2 F\g] with g := (w - 1) V™" 1 . 

Moreover, if r] := max{l - Wq, W\ — 1}, then lim ?? _ > o + (|1 — Pl\ + P2) = 0. 



The constant (3\ and p 2 are explicitly given in (15. 2p in terms of m, D*, Wq and W\. 



Proof. Define h k (w) := {w^ 1 - l)/(fc - 1). Let q := W, 
is non-decreasing, 



2(2-m) . 
j «1 •- 



Wl {2 ~ m) . Since \h 2 /h r 



(5.1) 



h' 2 (W ) 
h' m (Wo 



< 



h 2 (w) 
h m (w) 



< 



h 2 (W 1 ) 



< 



h 2 (W x ) 
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Note that oq = ao(Wo) < 1 < ol\ = ct\{W\) and both converge to 1 as Wo, W\ — ► 1. 
Using the fact that V^ 1-1 = + |x| 2 and an integration by part, we get 

[ \V [h k (w)V™- l ]\ 2 V D , dx 



\h' k (w)\' 2 \Vw\ 2 V^ 1 - 1 dx + 



1 — m 

m.2 



\x\ 2 \h k (w)\ 2 Vd, dx — d 



1 — m 



\h k {w)\*V% dx. 



Let g := (w — 1) 1 . Applied with k = 2 and k = m, the above identity gives 
![<?]= m / IV^HV^ 1 - 1 ]! 2 Vb. dx 



<mai / |Vw| 2 V; 



2 T /2m-l 



dx 



+ ai 



1 — m 



m J R d 



\x\ 2 \h m (w)\ 2 V Dt dx-d(l-m) / \h 2 (w)\ 2 V% dx 



and 



1 2 T A2m.-l 



dx 



= / |V [^(w)^- 1 ] | 2 V D , dx 

JR d 



1 — m 
m2 



\x\ 2 \h m (w)\ 2 Vb* dx + d 



1 — m 



IM«0r^D. dx 



Collecting these estimates, we obtain 
\[g]<m ai [ \V[h m (w)V^- l }\ 2 V D ,dx + d{l-m) [ { ai \h m (w)\ 2 - \h 2 (w)\ 2 ) V% dx . 

JR d JR d 



Note that 



, [ | V [h m {w)V r D n - 1 ] | 2 V D , dx < W Q ~ l J[w) . 

JR d 

Using F[g) = \ J Rd \g\ 2 V 2 ~ m dx = \ J Rd \h 2 (w)\ 2 V£ dx with g := (w - 1) U™"\ we obtain 

\[g]<PiJ[w}+p 2 F[w] , 

with 



(5.2) 



Pi ■-- 



ax _ W[ 



2(2-m) 



and [3 2 := 2 d (1 — m) 



W W 

Note that a = a (W ) < 1 < a\ = ai(Wi) and both tend to 1 as W , W\ — > 1. 



□ 



5.2. Entropy - entropy production inequality. 

Theorem 5.2 (Entropy - entropy production inequality). Let m £ (0,1), m ^ m*. For any 

function w = w{x) satisfying (H1")-(H2"), if 1 — Wo > and W\ — 1 > are small enough, then 
there exists a positive constant 7 such that 

lT[w] < J[w] . 

As we shall see in the proof, the constant 7 can be estimated as follows: 



m - C mtd d(l - m) 



(5.3) 



7 > 2 



( W± 

\W 



2(2-m) 



1 



C^dW^W, 



3 w 2(2-m) 
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The condition that 1 — Wq > and W\ — 1 > are small enough in the statement of Theorem 15. 2 
can be replaced by a weaker condition which amounts to ask that the right hand side of the above 
estimate is positive, that is, with X m ,d = m/Cm,d, 



i 



^1 < ( i + Xm < d 



W V d(l-m) 

Proof. Let g := (w - 1)V^ -1 . By definition of D*, = J Rd (w - l)V Dt dx = f Rd gV*- m dx if 
m > m*. By the spectral gap estimate (|4.6p (also see Theorem I A. ip . 

2F[g]<^f\[ 9 ]. 

By Lemma 15.11 

2F[g]<^\[g]<^L((3 1 J[ w ] + (3 2 F[g]) 1 
m m 

from which we deduce that 

F[g] < 9 Cm '^ x „ J"M • 

By Lemma the conclusion holds with 

2m- C m>d P 2 



7 



w^ l - 2 c m4 p 1 



under the condition 2m — C mj dP2 > 0. According to the definition f|5.2 j) of this amounts to 
the condition 

C m ,d/3 2 = 2(i(l-m)C m , d ( — -1) <2m, 

that is oii/otQ close enough to 1, which follows from the requirement that 1 — Wq > and 
W\ — 1 > are small enough. □ 

6. Convergence with rates 
6.1. Exponential decay of the relative entropy. 

Proposition 6.1. Let m E (0, 1), m ^ m*. There exists a positive constant 7 such that, for any 
solution w of (|2.2p with initial data u>o satisfying (H1")-(H2"), if 1 — Wo > and W\ — 1 > 
are small enough, then 

F[w(t)] < F[w ] e-T* . 



The value of 7 can be estimated from below by (J5J3J). 

Proof. We combine formula ^jF[ui(i)] = —J[w{t)\ with the estimate J[w(t)\ > ~f T[w{t)\ ob- 
tained in Theorem 15.21 and then integrate the resulting differential inequality. □ 

2d (1 — m) 

6.2. Moments, L p and C fc estimates. We recall that q* := — - and define 

a (1 — m) + 2 (2 — m) 

7(9) : = \ if Q G (<?*> 2 ] , 7(g) = f, + 1 N if g € ( 2 , 00) , l(oo):= 



2 " v " ' J ' g(d + 2) * v ' 7 ' ,v 7 d + 2 

The following lemma helps to understand better the consequences of the convergence of the free 
energy £ [v|VdJ> in terms of L p , moment and also C k convergence. 

Lemma 6.2. Let m G (0, 1) and consider a function v satisfying (H1')-(H2'). Then 
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(i) For any $ G [0, jEStL there exists a positive constant K$ such that 

x\*{v-V D ,)\ 2 <K* (£ [«|Vb.]) 1/2 • 

(ii) For any q G (q*,2], there exists a positive constant K{q) such that 

\\v-V D ,\\ q <K(q) (S[v\V D *}y {q) • 
Consider now a solution v of (11. 4ft such that vq satisfies (H1')-(H2') and fix some to > 0. Then 

(iii) For any j G N and any to > 0, there exists a positive constant Hj such that 

IK*) " Vb.|lcJ(R") < #j ( f K^bJ) 3 ^™ V t > t . 

(iv) For any q G (2, oo], i/iere exists a positive constant K{q) such that 

\\v(t)-V D J q <K(q) (£K*)|VbJ) 7(9) Vt>t - 
Proo/. (i) With := 2 su Pr . >0 r 2 * (D» + ^ r2 )-(2-m)/(i-m) j 

II |xr(« - Vb.) || 2 2 < J k* / |« - Vb. | 2 V™~ 2 dx < k$ F[(w - l)^- 1 ] , 

^ JR d 

and the right hand side is equivalent to £ [v\Vd,] by Lemma 13.11 

(ii) Let q G (g*,2). By Holder's inequality, 

/ \v-V D .\*dx=[ vt m)q/2 i\v-V D fV™- 2 ) q/2 dx<c(q)([ \v-V D fV™- 2 dxY 2 , 

jR d JR d * \JR d / 

(2— m) q 

where c(q) 2 ^ 2 ~ q ^ := J^Vp 2-9 drc is finite for any q > q*. By Lemma 13. 1\ the estimate 
holds with K(q) := c{q) l / q (2 W 1 2 " m ) 1 / 2 . In the limit case q = 2, the same method applies with 

c (2) = ||y D '- m) ||oo = D*- (2-m)/(1 ~ m) . 

(iii) We apply the interpolation inequality (|3.2|) to / = v(t) — Vb, and bound ||f(t) — Vb»||cj'(R' 1 ) 

in terms of |K£)||ctf+ 1 > which is uniformly bounded by Theorem 12.41 and \\v(t) — Vb, || 2 /^ +2 (- 7+1 ' \ 
for which we apply the result of Part (i) with # = 0. 

(iv) By Theorem l2.4l v(t) G Cl(R d ) and — Vd, is bounded in CI uniformly for any t > to > 0. 
By II, 

2d d 2 

|K<)-Vb.||oo < c^ 1/2 \\v(t)-v D j d c \i Rd) \\ v ( t )-v D Jt 2 ■ 

We conclude by using Holder's inequality, \\v(t) — Vb, \\ q < IKt) — Vb, ||oo ||^(i) — Vb, 1 1 2 - ^ 
6.3. Improvement of the convergence. 

Theorem 6.3. Let d > 3, m G (0, 1) tuii/i m 7^ m*. Consider a solution w of (12. 2ft u>it/i initial 
data satisfying (H1")-(H2"). There exist a positive constant K, and a time to > smc/i t/iat 

Moreover, for any A G (0, A mj d) ; t/iere exist a positive constant Coo & n d a time to > swc/i t/iat 

K(t) - l|oo < Coo e 2 2-™ Vt>t - 
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Here \ m ,d = m/Cm,d where C m ,d is given in Theorem I A. 11 Hence the rate of decay obtained 
by spectral methods for the linearized equation exactly gives the rate of decay for the nonlinear 
problem, and the price to be paid is only on the constant fC. As a subproduct of the proof, for 
some positive constants 770 and 700 which are defined below, we obtain the following estimate 



VQ (2-m) 

(t )] e t°°(i-"0 e 2A ™' dio . 



Proof. By Corollary 13.41 for any e > 0, there exists to > such that w(t) = w(t + to) satisfies 
Assumption (H2") at t = with < 1 — Wq < £ and < W\ — 1 < e. With e small enough, w 
enters in the framework of Proposition 16. 1\ with 7 as in Theorem 15.21 From now on, we assume 
that t > to and simply write w instead of w. 

On the one hand, by Lemma 16.21 and Proposition 16.11 we have 



\v(t) - VdJoo < o- e" 7oo{t " to) with a := if (00) F[wq\^ and 7 C 



7 



d + 2 ' 

which, in terms of it?, gives the estimate 



\w(t, x) — 1| < ctq e 



-7oo (t—to) 



1 

1 - m , , 2 



£>* + -z — \x 
2 m 



1 — m 

V t > t , V x eR d 



On the other hand, let h a {s) := (1 + s) a , a > 1. For any s G [0, so], 

h a (s)-l s u/lf , f« + f«(«-l)(l + ^or 2 if a>2, 

< a + — max h a (s) = < 

s 2 se [°' s °l [ a + f a(a-l) if a<2. 

Apply then this inequality with a = 1/(1 - m), s = (D* - Di)/(Di + |x| 2 ) < so = -§7 
to get the existence of a positive constant M.\ = Aii(m, D*, D\) such that 



1 

Similarly, for any s G [— sq,0] 



> a - — max h a {s) 

S 2 s6[-s ,0] 



a — 4^ a (a — 1) if a > 2 , 

a-fa(a-l)(l-s ) a - 2 if a<2, 



2 

so that, with a = 1/(1 - m), s = -(D - D*)/(D + |x| 2 ) < -s = ^ - 1, we get the 
existence of a positive constant .Mo = ■Mo(tti, D^^Dq) such that 



Hence, there exists a positive constant .A/f depending on max{.A/fo> A^i}, -Do an d D\, for which 
we obtain 

f an e -7oo (t-to) 1 i_ m , , 1 l=m 

\w(t,x) - 1| < min ^ ^ , M V^~ m \ = ae~^ 3=s (*-*») , CT : = Uq - m . 

As a consequence, for any t > to, we have improved bounds on w, with Wq and W\ replaced re- 
spectively by <7 (*) := l-cre" 700 ^ ( *" io) andcri(t) := l+o-e" 70053 ^ (i_io) . As in PropositionEH 
according to Theorem 15.21 and Inequality (|5.3p . z(t) = JT[tt;(t)] satisfies 

§ < -7(*)*(t) 
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with 



. ? n <7 (*) 2(2 - m) - C m . d d (1 - m) [a^tfV-™) - a (t) 2 ^] 



A TO ,d = rn/C m>d and r?(t) < r/ e 7 °° 2 
that for any t >to, 

, *(*) 

log 



for some t/q > 0. A Gronwall argument then shows 



'Jo 



i < — 2 A m d (t — t ) H r 



1-e 



< — 2 A m d (t — t ) + 



770 (2 - m) 
7oo (1 - m) 



which completes the estimate on F[w(t)). For t large enough, \~f(t) G (A, A m> d) and the L°° 
estimate follows. □ 



Proof of Theorem \l-4\ As in the proof of Corollary 13.51 a Holder interpolation inequality shows 
that, for any 5 > 0, 

/ \w(t) - l\ p dx < \\w(t) - if' 6 -** [ f-^-] 5+2 d x . 

jRd J R d \1 + \x\ 2 J 

□ 



6.4. Proof of Theorem 11.21 We first apply Theorem 16.31 and Lemma 16.21 (i), with i? = to 
obtain, for some to > 0, 

IK*) - Vb,|| 2 < K^a {£ [v(t)\V D ,]^ < C 2 e" A - d * V t > t , 

for some positive constant Ci- By the interpolation inequality (|3.ip with A = —ad < = fi < 
1/2 = u, C = C_ Q d, o,i/2; an d Lemma [3~6l (|3.4p . we have 

Kt)-Vb.||oo< -VbJ|&« ||u (t) - VbJ|£~ fl < cw 9 |Kt)-y D j|^ ||«(t) - VbJIa"* 

where = 1/(2 + ad). This implies 

IK*) - v D , Hoc < cva-*) % e /i i ~ e ) \\ v (t) - v »,h vt>t . 

From Holder's inequality, ||u(t) - Vb.||, < ||«(*) - VdJ|^" 2)/9 \\v(t) - V D J 2 2 ,q , q E (2,oo], we 
deduce that \\v(t) — Vb„|| g decays with the same rate as \\v(t) — VoJh- If Q £ (q*,2), we apply 
Lemma 16.21 (ii), and Theorem 16.31 to prove that for some positive constant C q and for some 
to > 0, 

Ht)-V D J q <C q e- x ^ dt Vt>t Q . 

Similarly, the estimate \\v(t) — VD,\\cj(R d ) follows from Lemma |6.2[ (iii), and Theorem 16. 3[ This 
completes the proof of Theorem 11.21 □ 



Appendix A: Hardy-Poincare inequalities 

In this appendix, we state and prove a result on inequalities which we have already been 
partially studied in [7]. Here we give more details and a few improvements. We are especially 
interested in the explicit values of the constants which enter in the convergence rates of Theo- 
rems 11.21 and 11.41 This is why we take weights which are adapted to Equation (II. 4|) and define 
the measures 

dfi := V^~ m dx and dv := Vd dx , 

where V D (x) = (D + ^ |x| 2 )" 1/(1_m) . Incidentally we observe that dfi = V^~ m du. To a 
function g £ Ll(d/i), we associate its average g = J Rd g(x) d/j,. Recall that = (d — 4)/(d — 2). 
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A.l. Statement and comments. 

Theorem A.l. Let d > 1 and D > 0. If m G (0, 1) and 1 < d < 4, or m G (m*, 1) and d > 5, 
then there exists a positive constant C m ^, which does not depend on D, such that 

(A.l) / \g-g\ 2 dfi<C m . d [ |V 5 | 2 aV V g £ V(R d ) , g=[ g dfi . 

JR d JR d JR d 

In case d > 5 and m G (0, m,), we /iawe 

(A.2) / g2dfi<C m>d [ \Vg\ 2 du V 9 GP(R d ) 

j 8m (1— m) . , • , 

and C m4 = [ (d _ 2 )( m - m .)j2 zs opizma/. 

Estimates of the optimal constant C m! d when m > m* are given below in Proposition IA.31 
With v m (x) = (1 + |rr| 2 ) _1 /( 1-m ), a simple change of variables shows that \ m ,d = ?n/Crn,d is such 
that 

J Rd | V/i| 2 Vd dx 1 - m / Rd | V/i| 2 v m dx 



(A. 3) A m fj = m inf — = ^ = inf — 

1 ' ' h f Rd \h-h\*V*- m dx 2 h yih-hfv^dx 

where the infima are taken over the set of smooth functions h such that 

- either m < m, and supp(/i) cR^ {0} and h = 0, h = 0, 

- or m > m*, 

- f Rd hV D m dx ^ ~ J Rd h v 2 m m dx 

Jr* m d% J^ d v m m dx 

This already shows that \ m ,d is independent of D. 

We observe that as |x| —> oo, dfi ~ dv/\x\ 2 . Hence, if m € (0, m»), Inequality (|A.2|) is of Hardy 
type. Otherwise, if m G (m*, 1), Inequality (lA.lj) involves an average and is rather of Poincare 
type. In such a case, we shall also say that it is a weighted Poincare inequality, or that there 
is a spectral gap, since for the associated operator, the lowest eigenvalue, 0, is achieved by the 
constant functions, and the second eigenvalue corresponds to \ m ,d = fn/C m ,d where C m< d is the 
best constant in the inequality. See [7J for further considerations on these issues. 

We also remark that Theorem provides an explicit example for which the weighted Poincare 
inequality holds, while the corresponding weighted logarithmic Sobolev inequality does not hold, 
even in dimension d = 1, as shown by [31 Theorem 3]. 

The proofs of (jA.ip and (|A.2p are quite different and for this reason we treat the two cases 
separately. We start with the proof of (|A.2p corresponding to the case m < m*, d > 5. 

A.2. Case m G (0, m*). The proof follows the ideas of [TJ. We reproduce it here for completeness. 
We compute 

\VV D (x)\ 2 = l -^-V D (x) 2 ^ 
ml 



and 

AV D (x 
V D {x) 

An integration by parts and the Cauchy-Schwarz inequality show that 



2m2 ~ ^ 3 _2m = 2dDm+(d-2) (m*-m)\x\ 



|ff| 2 AVb dx 



<2 / \g\\Vg\\VV D \dx 

JR d 

<2( f \g\ 2 \AV D \dx) f ( [ iV^lVVbHAVbl^dxV 

\JR d J \JR d J 
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As in [24], we remark that AVd has a constant sign and get the estimate 



\g\ 2 AV D dx 



[ \g\ 2 \AV D \dx < 4 f |V5| 2 |VVb| 2 |AVb| _1 da;. 

JR d JM. d 



Weights can be estimated on both sides of the inequality: 

|AVb| 2d Dm + (d - 2) (m* - m) \x\ 2 (d - 2) (m* - m) 



2-m 
D 

|VVb| 



< 



m(2Dm + (1 — m) \x\ 
2\x\ 2 



> 



< 



m (1 — m) 
2 



\AV D \V D ~ 2dDm+ (d - 2) (m* -m)\x\ 2 ~ {d - 2) (m* - m) ' 
which proves <fO]) . See [7] for further details. 

We now consider the limit D — > + . With a := l/(m — 1) £ (1 — d/2, —1), that is m G (0, m*), 
and 

8m(l-m) 1 — m_ 4(1 -m)2 
Kq "~ [(d-2) (m-m»)]2 ' 2m ~ [(d - 4) - (d - 2) m]2 ' 

Inequality (|A.2[) takes the form of a weighted Hardy inequality, 

2 



I \x\ a dx < K a f \Vg\ 2 \x\ a dx Vj£P( 

J9.d \X\ 2 jRd 



Such an inequality is easy to establish by the "completing the square method" as follows. Let 
a G R \ {q} with a* := 1 - d/2, and 5 G £>(]R d ). Then 







x\ 2a dx 



\x\ z 

\Vg\ 2 \x\ 2a dx+ X2 - \ (2a + d-2) 



fL k| 2a dx . 



An optimization of the right hand side with respect to A results in choosing A = (2 a + d — 2)/2, 
that is 

1 4 _ 4(l-m)2 

A2 ~ (d + 2a-2)2 ~ [(d - 4) - (d-2)m]2 ~ Ka ' 

The weighted Hardy inequality is optimal, with optimal constant K a , as follows by considering 
the test functions g £ (x) := min{e~ A , (|x| _A — e A )+} and letting e — > 0. □ 

A closer inspection of the proof reveals that the constant K a in the weighted Hardy inequality 
also is optimal when m > m*. Consider indeed the test functions g e (x) := \ x \ l ~ a ~ d l 2+e for 
\x\ < 1 and g e (x) = (2 — |x|) + for |x| > 1, and then let e — > 0. 

Proposition A. 2 (Weighted Hardy inequality). W%/t t/ie above notations, for any a G R, 



Ixl Q dx < K, 



R d \ x \ 



[ |V 5 | 2 |x| 2a dx Vg£V(R d ) 



with the additional requirement that g is supported in M. d \ {0} if a < a*, and K a is optimal. 

The range m G (0, 1) corresponds to l/(m — 1) = a G (—00, —1), so that m = m* is equivalent 
to ol = ca^ . Notice that the result holds without other restriction than but one has to be 

careful with integrability condition at x = if a < a*. 
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A. 3. Case max{0, 771*} < m < 1. Several partial results are known. In the range m £ (m c , 1), see 
[14] for an estimate of C m ^ based on the Bakry-Emery method, [7] for other estimates, and [27] 
for the exact values of the optimal constant for a corresponding linear problem. 

We now prove (fATl) with some explicit estimates of the constant C m d in the whole range 
(max{0, to*}, (d — 2)/(d — 1)) D (to* , m c ] . Because of the change of variables (|A.3|) . our task is 
now to characterize C m ^ as 

(1 -m)C m4 \~ 1 _ f Rd \Vh\ 2 v m dx 



1m J h J Rd \h-~h\ 2 v^ m dx ' 

On M + , consider the function fj,(r) := r d ~ l {\ + |r| 2 )( 2_m )/( m_1 ) , and denote its median by i]. Let 
v(r) := r d ~ l {l + r2) 1 /( m ~ 1 ) and define for all ( > the quantity 
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K(C) :=- max{A(C),B(C)} 

1 — to 



with 



A(C) := sup 



r 



n(s) ds 



? ds 



ds 



u(s) 



+oo 



fj,(s) ds 



B(C) := sup 

Jr v \ s )\ r>C UC, 

By convention, we take K(0) = ^ B(0). The following result is inspired by [3 [HJ QH [38] • 
Proposition A. 3. Let d > 1. For any m ^ m*, 

8m (1 — m) 
m ' d " [d - 4 - m (d - 2)] 2 ' 

If m £ (m*, 1), i/ien 

f 4 771 1 

C m ,i < K(0) and C m , d < max <^ 2 K(t/) , — — — } if d > 2 , 

I (1 - m) {d - 1) J 

where, for any m G (m*, (<i — 2)/(d — 1)), 

3 — 2??i / 2 — m 



77i (2 - m) 2 i-™ 1 + 2 i-" 1 

KfT)) < ^ 

Kl> - d [d-4-m(d-2)]2 

The function f^~ m is integrable for any m £ (m*, 1), so that K(t/) is well defined in this range. 
The upper bound on C m ^ is equal to its exact value up to a factor which is at least 1/4 (and at 
most 1). Such an interval is inherent to the method, see [38]. The case m = m c < (d— 2)/(d— 1) is 
covered, showing in particular that C mcj d is positive, finite. The bounds diverge as m \ m* with 
same behavior at first order. Our approach can be extended easily to the case ((d — 2)/(d — 1), 1), 
with slightly different estimates of K(t/), but this case is already covered in [T4"l [27] by other 
methods. The restriction m < (d — 2)/(d — 1) is convenient from a technical point of view, and 
not essential at all. 

Proof. The lower bound on C m ^ is achieved as in Section IA.2I by taking the limit D = 0, thus 
showing that C mt d > ^i/( m -i) an d using Proposition IA.21 

Let us prove the upper bounds. We introduce the standard change of variables from Cartesian 
to spherical coordinates, i.e. r = |x|, and •& = x/\x\. In these coordinates, the gradient can be 
written as (d r , jVg) where d r = ~ ■ V is the partial derivative with respect to the radial variable 
r and Vg is the derivative with respect to the angular variables. We shall denote by S^ 1 C M. d 
the unit sphere and parametrize it with the variable "d. 

The radial density functions r i— > /i(r) and r i— ► u{r) are such that v m dx = [x{\x\)dx and 
v^ m dx = dx. We introduce the following normalization constants: 

r 2ir d / 2 — r 



26 A. BLANCHET, M. BONFORTE, J. DOLBEAULT, G. GRILLO, AND J.L. VAZQUEZ 

With these notations, 



fj,(r) dr d$ 



v m {x) 



1 + \x\ 2 

We define a directional average of a function / by 



+00 



dx = v m 771 dx and v(r) dr dt? = v m dx 

H(r) 



/(r, $) fi{r) dr with fi(r) :- 



Jo 00 K s ) ds 



and the global average of / by 

JrcI f m dx 



f:-- 



2-m 

Jm 



dx 



f{r, 1?) /t(r) dr dt? 

(0,oo) xS 1 *" 1 



In the case d = 1, Theorem 2 of [1], also see [38], says that Inequality (jA.ip holds with 



2 m 

C m ,i < sup 

1 - m r>0 



+00 



/i(r) dr 



dr 



v r 



K(0) 



in which case we also have the estimate C m> i < K(0) < 4C m> i. 

In case of radial functions, Inequality (jA.ip takes the form: 
+00 _ 2 j _ m 



(A.4) 
with 



/(r)2-/ n(r)dr< 



2 m 



C* d < - — max <; sup 



1 — m 



r>rj 



fi(r) dr 



2 m 

dr 

v(r) 



/^rad 



sup 

r<rj 



|/'(r)| u{r) dr 
"> dr 



n{r) dr 



. v{r) 



K(r?). 



It is straightforward to show that K(7y) is finite, with the present choices of \i and v, for m £ 
(m*, (d - 2)/(d - 1)), and as above, C™ d d < K(rj) < 4C™ d d . 

We now focus on the case of non radial functions, with d > 2, and rewrite the left hand side 
of (jA.ip in spherical coordinates. 

2 - rr . ~ |2 



/(*) - / 



2-m 
- 01 



dx 



= 

< 2w d [(I) + (II)] 



|/(r,#) - / I //(r)drdtf 

(0,oo)xS d " 1 

|/(r, ,?) - ^(0) + ^(0) - J| 2 M (r) dr &d 

(O^JxSd- 1 



with 



(I) 
(II) 



|/(r,tf)-/ M (tf)|V(r)drdtf , 

(0,oo)x§ d - 1 



/^(r)drd?? 

(0,oo)xS d " 1 



We estimate (I) by (|A.4p and get 



(I) 



|/(r)2-/ M (i?)|^(r)drd^< 



(0,oo)xS d - 1 



1 — m 
^^— 

2 m 



->rad 

■"m,d 



|<9 r /(r, 0)| z/(r)drd$. 

(0,oo)xS d ~ 1 



To estimate (II), we rely on the Poincare inequality on the unit sphere § 



d-l 



\u — u I 2 d$ < 



d-l 



\Vi)u\ 2 d$ VueHl(§ 



d-1^ 



Here u := Ld-i ud$. In the inequality, l/(d — 1) is the optimal constant, as can be checked using 
spherical harmonic functions. See for instance [5j El [43] . The inequality itself can be recovered by 



ASYMPTOTICS OF THE FAST DIFFUSION EQUATION 



27 



various methods. For example, using the inverse stereographic projection, see [37], the optimal 
Sobolev inequality on IR'' becomes 



\v\ p dd\ < 



gd-i a — 1 



\V#v\ 2 d# , 



for any u £ with p = 2d/(d — 2), d > 3. The inequality also holds true for any 

p G (2, 2d/(d — 2)) if d > 3 and for any p > 2 if d = 2, see (BJ. Hence we recover the Poincare 
inequality on S d_1 by writing v = 1 + eu and keeping only the terms of order e2 as e — > 0. 



We apply the Poincare inequality with u = 



fu f 



d- 1 



dtf 



Recall that |V/| 2 = |5 r /(r, #)| 2 + ^ |V^ /(r, #)| 2 . Using the Cauchy-Schwarz inequality and the 
estimate r2 /i(r) dr < v(r) dr, we get 



dl? 



< 



|V # /(r, ti)\ n(r)dr&d 



(0,oo)xS< 1 - 1 



1 |Vtf f(r,$)\ 2 i/(r)drch9 . 



r2 



This proves that 



(II) < 



1 



(0,oo)xS 



— |Vtf /(r,i9)| 2 i/(r)drdT? . 
d-i r2 



d- 1 

Summarizing, we have shown that 
By undoing the change of coordinates as in (|A.3[) . we get 



1 



|/(x)-/rd/x<max<!2C d d 



4m 



|V/| 2 7; m dx. 



|V/| 2 d^ 



(1 - m) (d - 1) 

The bounds on K(ry) follow by quite long but straightforward calculations, omitted here. 



□ 



Appendix B: Extension to exponents m < 

The presence at several instances of factors of the form 1/m in the previous calculations may 
suggest that there is an essential divergence as m — > 0. In the present section we want to dispel 
that impression by introducing a normalization that is often used in the literature, consisting in 
a rescaling of the time variable of the form r' = mr that modifies equation (jl.ip into 

(B.l) d T ,u = V • (u m_1 Vtt) ■ 

One of the first consequences is that the new equation, that we will call modified fast diffusion 
equation for clarity, makes perfect sense as a nonlinear parabolic equation of singular type for 
all the range of exponents m € R (including m = 0), in particular for all m < 1 that form the 
extended range of the fast diffusion. Such approach has been consistently used in [U] where it 
is shown that the effect on the self-similar solutions of Barenblatt type is just to eliminate the 
denominator m in the formulas (jl.2p . (|1.5p . Note the rescaling to obtain (|B.ip from the standard 
fast diffusion equation (jl.ip when m > can also be done by changing the space variable in the 
form x = \/mx' and not changing time. 

Since the general theory (existence, uniqueness, estimates, special solutions and extinction) 
has been developed to the measure we need it, we can follow the different stages of the present 
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paper with due attention to chasing the m factors, and the results stated in Section [T] remain 
valid. For instance, the formula defining J-[w] in Section 12.51 has to be replaced by 



In the linearization of Section |4] there is no m factor in the definition of operator A m and neither 
in the definition of \[g]. Let us mention two other points of interest: the exponent m* becomes 
negative for d = 1,2,3 and zero for d = 4, but it still plays the same role of an important critical 
exponent separating different behavior types. On the other hand, the constant \ m ,d that gives the 
decay rate in our main result has a finite positive value as m — > 0, according to Theorem lA.il This 
constant determines the rates in all results concerning the asymptotic behavior of the solutions 
and is not affected by our m-rescaling. 

We have refrained from treating the extension to m < in the main body of the paper in 
order to avoid further distractions in an already very technical matter. Whole details will appear 
separately. 
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